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Abstract 



^ ' The response of isovector and isoscalar pairing to generalized rotation in isospace is studied. 

OO ■ Analytical expressions for non-selfconsistent solutions for different limiting cases of the model 

o 

are derived. In particular, the connections between gauge relations among pairing gaps and the 
' position of the isocranking axis are investigated in N=Z nuclei. The two domains of collective 

and non-collective rotation in space are generalized to isospace. The amplitudes for pair-transfer 
^ , of T = and T = 1 pairs are also calculated. It is shown that the structure of the T = state in 

^ I odd-odd nuclei prevents any enhancement of pairing transfer also in the presence of strong T = 



correlations. The energy differences of the T = and T = 1 excitations in odd-odd nuclei are 
qualitatively reproduced by Total Routhian Surface calculations. 
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I. INTRODUCTION 



The event of detailed experimental studies of heavy nuclei along the N=Z line has resulted 
in a revival of theoretical investigations related to the properties of proton-neutron {pn) 
pairing correlations |]I|, ^ ||, ^, |], 0, |, |10|, |Tl|, |12|. In N=Z nuclei one expects short 
range correlations of pn type to be of importance due to the four fold degeneracy of the 
nuclear wave functions. Of particular interest is the role played by t=0 pairing correlations. 
It is still an issue of debate to what extent the correlations in this channel of the pairing 
force are characterized by a static gap similar to that of the well established t=l seniority 
force. 

The aim of our work is to devise a mean field model in which we can describe consistently 
excitations in real- and iso-space on the same footing. In a serie of papers we have shown 
that the isobaric analogue states as well as ground state masses along the N=Z line form 



an unique probe to t=0 pairing correlations p, 0, |T^, |T^. Within the mean field 
approach, the energy of isobaric analogue states can be described by means of the isocranking 
approximation, analoguous to the rotational excitations in real space [1^, [1^ . The associated 
broken symmetry is the deformation of the pairing field. The direction of the pairing field 
vector is in turn intimately linked to the direction of the cranking axis in isospace. Hereafter 
(Sec. fTT| and [TV] ) we derive analytical expressions to elucidate these relations which, in 
fact, determine the regimes of collective and non-collective rotations in isospace. These 
expressions and conclusions may be of potential use in other, yet unknown, double-phase 
paired systems that can be described by means of an external cranking-type hamiltonian. 



The collectivity of pairing correlations can be acessed by means of pair transfer In 
the analysis of nuclei in the vicinity of the N=Z line, it was concluded that indeed, the 
t=l, pri-pairing exhibits collectivity [|l^. From this analysis, it was suggested that there is 
little evidence of t=0 collectivity 1201 . Let us point out that similar conclusions concerning 



the role of t=0 pairing has been drawn by Berkeley group ||2T| based on the analysis of 
excitation energy spectra in N=Z nuclei. However, to draw definite conclusions a detailed 
understanding of the structure of both T = and T = 1 ground states in the odd-odd 
(o-o) nuclei as well as the structure even-even vacuum are necessary. There are empirical 
arguments based on isobaric symmetry as well as theoretical arguments based essentially 
on time-reversal symmetry breaking, which indicate that the structure of the T = 1 and 
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T = states in o-o N=Z nuclei is entirely different |T^. The simplest scenario, within the 
BCS theory, consistent with the data can be reached by assuming that the wave function 
of the T = state in o-o nuclei is a two quasi-particle excitation (2QP), whereas the T = 
ground state in even-even (e-e) nuclei and the lowest T = 1 state in o-o nuclei are both 
local quasiparticle vaccua i.e. OQP states |T^. Within this interpretation the transfer of 
a T = deuteron-like pair will always be strongly quenched, irrespectively of the strength 
of the T = correlations, see Sec. Hence, pair transfer may not necesserily be a good 
indicator for the strength of T = pairing correlations. 



The outline of the paper is the following: In Sec. |T| we derive analytical solutions for 



the isocranked model with t=l pairing only. In Sec. |^ we extend these solutions to include 
also t=0 pairing. In Sec. we discuss the role of the isospin symmetry breaking mechanism 



due to number projection and its influence on the Wigner energy. In Sec. VI we discuss pair 



transfer from T = and T = 1 states in N=Z o-o nuclei. In Sec. [VI I| we investigate the 
influence t=0 pairing on nuclear deformation by performing Total Routhian Surface (TRS) 
type calcualtions for T = and T = 1 states in N=Z o-o nuclei. In Sec. |VIII| we discuss 
briefly the shortcomings of our model due to the lack of the particle-hole isovector field. We 
summarize and conclude the paper in Sec. ^X. 



II. GAUGE INVARIANCE PROPERTIES OF PROTON-NEUTRON COUPLED 
HFB EQUATIONS 

(")=(")• 

V — > e''^V (2) 

which we call later global gauge invariance transformation (GGI) since it requires both 
proton Vp e'^^Vp and neutron Vn — > e'^^Vn amplitudes to be transformed simoultanously. 
Indeed, in this case the density matrix and pairing tensor transform as: 

p = V*V^ — > p and K = V*U^ — > e^^'^K. (3) 
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The HFB (BCS) equations: 




are invariant under the transformation: 



Hence, the single-particle potential F oc p and pairing potential A oc k; become 

r — ^ r and A — ^ e'^'^A. (4) 

The GGI allows, for example, to choose one of the pairing gaps to be real and we will take 
advantage of it by assuming [apart of Sec. |rV|] that the neutron gap A„„=A>0. However, 
in some cases and for sake of simplicity, the eigenvectors will be given only up to the gauge 
transformation. 



III. TWO-DIMENSIONAL ISO-CRANKING SOLUTIONS OF THE t = 1 PAIR- 
ING MODEL 

Let us consider in this section a model hamiltonian [h=l for convenience] : 

H^^ = hsp - Gt=iPlPi - (^X (5) 

containing an iso-symmetric particle- hole (ph) mean-field hgp. The isovector t=l pairing 
interaction is generated by: 

A±i = 5^ aL(p)C(p) and Pl^ = ^J2ialal + alal), (6) 

and we consider two-dimensional iso-rotations a3T-=[u;^ cos cu^ siny), 0]. Planar rotation in 
N=Z nuclei are the most general since (f^) = due to number conservation. We are 
interested in analytical but non-selfconsistent BCS solutions within the constant gap ap- 
proximation. Moreover, we assume that A„ = Ap = A and |A„„| = A = |App| but do not 
make any further restrictions concerning the phase relations between the gaps. Within the 
BCS approximation and under the above mentioned constraints the problem reduces to a 
diagonalisation of a 4 x 4 matrix with A and A^ being real [cj = Cj — A] : 

' ~ei-Ei -\ujrt-'^ A e^'^Ao 1 I" f/^,„ 
_ e-'^Ao e-^"A ^u^e-'^ -e, - E^ \ [ Vi,p 



(7) 
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The physical (positive) eigenvalues of are: 

Ei± = 

where 



A2 + A2 ± VXi 



(8) 



X, = e^col + AA'Al cos' - f ) + co^A' sin^ - |) 

-4eia;^AAo cos (^ip - ^ j cos (^^/^ - (9) 

These roots are double-degenerated [Kramers degeneracy] with eigenvectors of the form: 














Ui±,p 






















U~i±^p 











Vi±,n 













































(10) 



A. Model including t=l, ±^=^1 pairing 

It is pedagogical to consider solutions to the model step by step starting with the standard 
case of a pn unpaired system Aq = 0. In this case the eigenvalues reduce to: 



Ei± = J A^ cos^ (^(f 



UJr 

y 



±Jej + A? sin^ [ip - 



(11) 



Searching for the eigenvectors we assume that: (z) Satisfactory solutions must obey a min- 
imal consistency condition and reproduce the initial [see eq. (0)] relative phase relation for 
the gaps: A„„ = e*"App. (m) Moreover, we assume that neutron and proton amplitudes can 
differ only by their phase i.e. 



e'^Ui^p and Vi^n = e'Wi^p 



(12) 



which appears to be a reasonable assumption especially for systems described by common 
state independent order parameter A. Under these conditions we encounter only two types 
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of solutions: 



^-- = {2k + l)- 
(f-- = {2k)- 



(13) 
(14) 



2 ' '2 

In the first case ip — a/2 = {2k + l)vr/2 the quasi-particle (qp) routhians ([TTD become 
hnear as a function of iso-frequency: 

1 



Ej ± -u^ 
2 



where Ej = we| + A^. 



(15) 



Before the first band crossing i.e. for frequencies lower than 

Ur < 2Ei = 

where Ei denotes energy of the lowest qp at frequency zero, the qp routhians are 

Ej± = Ej ± -Ur. 

The associated eigenvectors are also independent on Ur and equal: 



(16) 



(17) 



Uj±,n 










Te'^Uj 












_ ±e-^^Vj _ 



(18) 



where the Vj and Uj amplitudes are: 



For frequencies cJr^'' < ujr < ujY'' = 2E2 etc. standard procedure for blocked states can be 



,(2) 



_ 1 

'~2 



A 



(19) 



applied to calculate eigenvectors, see e.g. [22 



The situation described here is characteristic for non- collective type of rotation. Indeed, 
in this case the iso-alignment will change stepwise at each crossing frequency similar to 
the sp-model described in detail in Refs. W\ ITR |TB|. 



The second possibility (|T4D ip — Oi/2 = k-n yields routhians of standard BCS-type: 



(20) 



where ej± = ej ± ^Ur- The eigenvectors are: 



(21) 



where 



1- 



and Uj± 



1 



(22) 



In this case the qp-routhians have a nontrivial dependence on lOr- They give rise to a 
smooth ahgnment processes typical for collective rotation. The particular case of this class 
of solutions corresponding to ^ 



a 



was discussed in detail in ref. |jT6 
These two classes of solutions have a simple geometrical interpretation+cite[Gin68,Fra99] 
Let us define the vector of anisotropy in isospace A = [A^., A^, A^] as: 

1 



In our case Aq = and: 



— I 

7i 



(App + A^„„); A, = A, 



pn 



A. 



— tan 



a 



(23) 



(24) 



as shown schematically in fig. |l]. Relations (|T3D-(|14D position the axis of iso-rotation either 
parallel or perpendicular with repect to A giving rise to non-collective or collective iso- 
rotation, respectively. In particular, for most standard choices of phases A„„ = ±App the 
collective axis is the x(?/)-axis, respectively. One may summarize that the possible solutions 
in this case only allow for principal axis cranking, where the cranking axis is determined by 
the gauge angle of the neutron and proton pair gap, respectively. 



B. Model including complete t=l pairing 

According to the geometrical interpretation of the pairing gaps, switching on adiabatically 
pn pairing Ap„[= A^] should always induce collectivity. Indeed, in this case the qp routhians 
(§) would depend on uj-r in a complicated, nonlinear way for both cases (|1^) and (|14|). Closer 
examination shows, however, that the non- collective solution (|THp cannot be generalized 
to accomodate pn pairing. On the contrary, collective solutions (^Tj) can be rather 
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FIG. 1: Schematic drawing showing the relative position of the axis of anisotropy [A] versus 
axis of isorotation. The phase relations (|l3|)-(|l^) allow either for parallel [non-collective case] or 
perpendicular [collective case] position of the rotation axes only. Note that tilted solutions are 
not allowed and hence, three-dimensional isocranking in N=Z nuclei can be effectively reduced to 
one-dimensional theory. 

strightforwardly extended to include pn pairing. In this case the qp routhians take the form: 

(25) 



E. 



I A =F Aoe*('/'-'^)|2. 



and the eigenvectors are: 



(26) 



where 



E. 



and Uj± = - 



1 A T Ap„e-*'^ 



(27) 



Let us observe that for — ip = [2n + l)7r/2 the routhians (|25|) have a similar form like 



discussed by Goodman ||2^ for the static case of ^,-=0: 

Ei± = 



ej^ + A2 + A2. 



(28) 
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For ip — (f = nvr, on the other hand, one gets routhians similar in structure to those discussed 



by Bes et al. |20 



E^± = + (A T (-1)"A„)^ (29) 

To gain further inside into our solutions let us calculate the density matrix: 

Pir,iT = Pir^iT = -A'^ - ^ - ^\ for T = n,p (30) 













ei- 1 







4 

Pm,ip = Pin'ip = ~, \ ~^ ^ f" (31) 



and pairing tensor: 



l^in,ip = l^ip,in = 7^0^ '^^j ^ ~ "^^j ^ (33) 



4 

where 



xW = -i-±-i- and XW=5:X(^) (34) 



Then, the gap equations are: 



^ - ^e'(^-^)x(-) (35) 



4 - Ae-(^-^)x(-). (36) 



G't=i,i2=o Aq 

For isospin symmetric pairing Gt=i,t:,=±i = Gt=i,t^=o = G'f=i the equations can be solved 
either when: (a) X*-"-* = and for essentially arbitrary values of gaps or (b) for A = Aq 
and %l) — if = niT. The latter case leads to gapless superconductivity ( ^9|) with a rather 
unphysical qp-spectrum. 

In particular, for the equidistant = i6e level model with symmetric cut-off, the particle- 
hole symmetry leads to X^") = but only for the routhians of Goodman- type (p8|). In such 
a case the particle number condition X = Trp is also automatically satisfied since [see eq. 

(H)|: 

The alignment in iso-space is then given by the formula: 

T^/&i^T^4E{|;-H (38) 

j>0 ^ ^ ■' 



which is analogous to the one obtained in ref. [there, should be formally replaced by 
A^ + A^] . Moreover, since for the isospin symmetric model A^ + A^ is constant the moment 
of inertia in isospace (Mol-i) and all conclusions drawn there remain unaffected. 

For more complex routhians (e.g. given by eq. (^) we have not been able find analytical 
solutions even for such high-symmetry like the equidistant level model. 



IV. MODELS INCLUDING ISOSCALAR PAIRING. 
A. Pure t=0 pairing model. 

Let us now consider the model hamiltonian: 

H'^^ = hsp — Gt=iPiPi — G't=o-fcl-Fo ~ '^rt. 



(39) 



consisting of both isovector and isoscalar pairing interactions coupling particles in time 
reversed orbits: 



>t 



V2 



(40) 



i>0 



Let us start our considerations with a pure isoscalar t = model. In this case the BCS 
equations take the following form [the pairing gap Aq can be chosen real]: 



Ci - Ei 






-Ao 




ei - Ei 


Ao 






A* 


-ei - Ei 


\ujre^^ 


—A* 

^0 




\uJre~'^ 


-h - Ei _ 



i,n 



Vi 



i,n 



(41) 



The eigenvalues are linear in uj^ and equal: 



E,. 



i± 



± Ei 



Ei± -ujr where Ei 



~4 + Ag. 



(42) 
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The equations are valid below the first crossing frequency uOr < ^'^^ = 2-Ei. The roots are 
double-degenerated due to the Kramers degeneracy and equal: 



where 












u. 


Ui±^p 











Ui±^n 






















'^i±,n 




V, 



























V. 

















and Vi 


2\ E, 



(43) 



(44) 

independently on Ur- The situation is similar to that described by (|T5|)-(|l8|) for the case of 
t=l, £^=±1 pairing in Sec. ill A . However, since the pair field is isotropic in isospace only 
non-collective iso-rotation takes place regardless of the direction of the iso-cranking axis. 



B. The t=0 plus t=l, tz=±l pairing model. 

The extension of the pure t=0 model including t=l, tj,=±l can be done by either linking 
it smoothly to non- collective solutions (|18|) or collective solutions (0). Let us consider first 
non- collective solutions obeying the condition (|13]). These can be relatively easy generalized 
to include t=0 pairing provided that the following phase relation is satisfied: 

9 -^a = rni (45) 

where Aq = e*^|Ao|. This is a necessary condition for the linear term of the determinant of 
the BCS matrix to vanish and, in turn, to preserve the mirror symmetric structure [Ei, —Ei] 
of the HFB solutions. In this case our solutions have non- collective character. The quasi- 
particle routhians below the first crossing frequency take the following form: 

E,± = E,± where Ei = e] + + \Ao\^ (46) 
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The associated eigenvectors are: 





























Ui± 


















^±,n 































Vi± 














where 



1 ATe-'^Ao 



However, since the pairing tensors are: 



and K. 



Ao 



(47) 



(48) 



(49) 



2^ ■ -m,ip ■-vp,m 2E- 

the selfconsitency conditions for pairing gaps yield the additional constraint: 

for the coupling constants. Solutions are therefore possible only for Gt=i=Gt=o- Coexistance 
of this type was already reported in the literature [|I], 0, ^. 



Generalization of the collective solution ([131) to include t=0 pairing is far more difficult. 
It can be shown, however, that for the following phase relations: 



the eigenvectors take the following form: 

























Ui±^n 











Ui±,p 











Vi±,n 




Vi± 









































(51) 



(52) 
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In this case the quasiparticle routhians are: 

E.± = (^e1 + ^1 + + |AoP ± 2^^(e1 + |Ao|2)+A2|Ao|2j (53) 

For the general case of cj^ 7^ the amphtudes are given by lenghty and rather non- 
transparent expressions. However, the property of these solutions are easily recognized 
already at Ur = when: 



e 2 / Cj 09 — 6' . ^, 

Ui± = 1= , COS (55) 

72 VE,±(E,±-e.) 2 ^ ^ 



and 



^^± = V^? + (^±l^o|)^. (56) 
The gap equations take the following form: 

4 ^(^) ^ ^x(-) (57) 



^^=1 A 
7^ = ^^"^ + ^^^-^ (58) 

<J't=0 1^0 1 

where X*^^-* are defined as in eq. (0). Since X*^^) < the gap equations can be solved only 
when Gt=Q=Gt=i and A = |Ao|. In this case the quasiparticle spectrum takes the rather 
unphysical gapless form in spite of the strength of pairing force. Let us mention here that 
particle-number is in this case automatically satisfied at least for equidistant spectrum since 
in this case particle-hole symmetry [eq. (|37|) ] is fulfilled. 

V. ISOSPIN SYMMETRY BREAKING 



The considerations of Sec. |ITI| show that collective iso-rotations are possible when the 
isovector pairing field A is perpendicular to the axis of iso-cranking. Such a collective 
motion allows to take into account isospin fluctuations in a static way by replacing the 
cranking constraint: 



;(t.)^ + (g^ = T (59) 

by 

(i)^ + {iy)^ = VT{T+l) (60) 
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In this way the Wigner energy [Ew ~ T] can be, at least partially, restored. Quantitative 



calculations [15], 0] show, however, that in the presence of the standard t=l field, 
the Mol is too large and cannot account for the empirical data. In these calculations the 
isovector part of the particle-hole field is not taken into account (see also Sec. |V111D . In 
the absence of these correlations the mechanism which lowers the Mol in isospace was fully 
ascribed to isoscalar pairing in complete analogy to the mechanism of lowering the spatial 
Mol by the isovector superfiuidity. 



However, as demonstrated in Sect. |^ at the level of the BCS approximation one cannot 



obtain a mixed collective solution which in turn is necessary to lower the Mol and to account 



for isospin fiuctuations in a static way (|60|) . To make use of the static condition (^) and 
to restore (locally) the Mol in isospace it is necessary to find phase mixed solutions. It has 
already been shown [Q] that approximate number projection of Lipkin-Nogami type allows 
for the mixing of t=l and t=0 pairing phases, breaking the isospin symmetry. The mixing 
appears only when the strength Gt=o of the average isoscalar field exceeds the strength of 
the isovector pair field Gt=i i.e. for = Gt=o/Gt=i > Xcrit ~ 1-1 |3 '^^^ Wigner energy 

(2) (2) 

generated through this mechanism contributes with opposite sign to Xtt and A^_^ building 
up an asymmetry in the auxilliary LN fields: 

X?iANrANr,. (61) 

Since the LN solution does not allow for an isovector pra-field, the situation qualitatively 
resembles the case considered in Sec. [IV B| . Let us however point out that the LN procedure 



introduces an explicit state dependence of the effective gaps 

^2) 

Moreover, the Lipkin Nogami parameters A^^, are anisotropic in isospace. In turn, a wider 
class of solutions become possible, including solutions which mix t=l t=0 pairing phases 
and are collective in isospace. The absence of the isovector pn-field in the LN solutions 
limits, however, the possibilities to calculate the amplitude of the isovector pra-pair transfer, 
see next section. 
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false vacuum 



even-even 



odd-odd 



A+2 



FIG. 2: Schematic figure of the pair transfer and structure of the T=0 and T=l states in o-o 
{A + 2) and ground state of e-e {A) N=Z nuclei. The thin arrows indicate the structure of the 
T=0(1) states in o-o nuclei in our model. Note, that the different structures of the T=0 states 
imply the quenching of isoscalar pair transfer even for a t=0 paired systems. 

VI. PROTON-NEUTRON PAIR TRANSFER 



The theory of pair transfer/stripping processes like (a,d),(^He,n) etc. were developed 



already in the end of 50's and the beginning of 60's |2^, ^ based on the plane wave 
Born approximation. The general expression for two-nucleon spectroscopic factor (2nSF) 
carrying nuclear structure information was derived using the spherical shell model in the 
jj— coupling limit, see Refs. p7|, |28|. The first analysis of the 2nSF using the pairing in- 



teraction model was done by Yoshida |2^. He pointed out the possibility of an enhanced 
cross-section for the two-particle transfer of the isovector pair due to collective pairing phe- 
nomena which may be particularly strong if the pairing coherency extends over few j— shells. 
Based on Yoshida's work, Frobrich ^ has analysed the influence of pn pairing on pn- 
pair transfer in N=Z nuclei using both t=l and t=0 pairing interactions within a single- j' 
shell model space. He pointed out that pn-pairing can enhance the cross-section by a factor 
of 3 as compared to conventional shell-model calculations of [^ . 

In the deformed shell model [paired mean-field] the differential cross-section describing 
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the transfer of a structureless isoscalar or isovector pn pair is proportional to: 



da 

^ ~ 1^10(00) I - l\^A|^10(00) 



^loloor^P = I(*a|Pio(oo)I^a+2(T))P (63) 



where we have assumed that the transfer goes from the T = 1(0) ground state \^ a+2{T)) 
of the o-o N=Z={A + 2)/2 nucleus to the ground state of the even-even nucleus 



N=Z=A/2. Within our model [|T4|, [T^, |T6| I'^a) has the structure of a Oqp state while 



the structure of |\E'^+2(T)) depends on the isospin T. For T=0, \'^a+2(T = 0)) is a two- 
quasiparticle state while the \'^a+2(T = 1)) states maintain the Oqp structure but is cranked 
in isospace, see fig. ^ In all cases the I^^a) and \'^a+2(T)) states are described by the 
fully self-consistent amplitudes (L''*-'^'' V^^^) and ([/(^+2) Therefore, one can make 

use of the generalized Wick's theorem to derive the explicit expression for the pair transfer 
amplitude T(^'^+2) (f 



Tl:^^ = ^(4.|4.«> J: ± (64) 



i>0 

where 



The overlap is given by Onishi formula: 



(^a|^a+2) = VBetU. (67) 

According to our model there is a fundamental difference between the structure of the 
T=0 and T=l states in the parent [o-o] nucleus, see fig. 0. The T=0 states correspond to 
2QP excitations. Therefore the deuteron will be transferred from a 2QP state to the BCS 
vacuum in the daughter nucleus. In such a case the Tq^'^^'^^ amplitude will essentially probe 
~ {Ua)^ [i enumerates blocked qp state] and will always be severly quenchend . 

On the other hand the isovector pn-pair will be transfered between Oqp states and can, 
in principle, be enhanced due to the pairing collectivity. The generic situation is illustrated 
in fig. ^ for the case of ^^Cr and ^''Mn. In these calculations we assume = 1.3 for both 
the ground state of ^^Cr and the false vacuum of ^°Mn. With increasing Ur ^°Mn picks up a 
small fraction of the isovector pn pairing and T^q'^'^'^^ increases reaching maximum around 
Ta; ~ 1. Then in the region of the phase transition it decreases rapidly due to the vanishing 
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l<|Pio|T=l>| 




0.5 1.0 1.5 2.0 2.5 



FIG. 3: Isovector pn-pair transfer {\Piq\uJt-) [upper part] and isoalignment [lower part] versus ujr- 
The calculations were done in the LN approximation. In this approximation t=l pn pairing plays 



a redundant role [0, |13|] resulting in the quenching of (|Pio|i^t-)- 



overlap, see eq. (Q). In the region of ~1 to ~ y/2 the |T|*-|^''^'''^-*| ~ l^ljo^ '^'''^''l ~ 0.4 i.e. 
no enhancement is calculated for the pair transfer. 

The quenching of the T^q''^^'^^ amplitude is related to the fact that within the LN model 
the isovector pn pairing plays essentially a redundant role. To investigate the importance 
of isovector pn correlations we have to come back to the t=l pairing model of Sec. UIB . In 
this case an analytical expression for the transfer amplitude can be easily derived: 

-,{A,A+2) 
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i>0 



(A)* (A+2) 



(A)* (A+2) 

u]_' vL 



X, 



(-) 



(+) 



where 



X 



(±) 



(A)* (A+2) (A)* (A+2) 



(68) 



(69) 
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The amplitudes {u, v) define the Bogolyubov transformation in the canonical basis: 

\±) = ^{\n)±e^\p)) (70) 

^) = -^(^)±e^>)). (71) 

They are equal: 

Ui± = V2U*^ and Vi± = V2V*^ (72) 

where Ui±, Vi± are given by eq. (^7|) [for the daughter nucleus they are obtained at uJr=Q]. 
Let us now assume that ^7^=0. When Ap„=0 then both m+ = n_ and f+ = f „ and we have 



an exact cancellation of two relatively large terms in the r.h.s of eq. (|68|) . On the other hand, 
for Ap„ 7^0, also ^ U- and we have a positive interference originating entirely from the 
terms proportional to Ap„ and, in turn, a rapid increase of IT^"^'"^"*"^)! as a function of the 
contribution of the pn pair gap to the total gap A|. = A^ + |Ap„p = const. This is shown 
in fig. ^ in version a) In this version of the calculations [|Ap„p = a^^A^] was forced to 
be the same for ^^Cr and ^°Mn (at Ur = 0). Version b), on the other hand, assumes fixed 
structure A = |Ap„| in ^^Cr. In this case refers to ^''Mn. It is clearly seen from the 
figure that the maximum transfer is expected for similar [but no necesserly equal] content 
of pn pairing in the parent and daughter nucleus. Obviously, other factors like different 
deformations etc. may additionally hinder pra-transfer. 

VII. TOTAL ROUTHIAN SURFACE CALCULATIONS INCLUDING 
ISOSCALAR PAIRING 

We have performed systematic calculations of the energy differences AEt = Et=o — Et=i 
between the lowest T = and T = 1 states in o-o N=Z nuclei using the Total Routhian 
Surface (TRS) method. We used the deformed Woods-Saxon potential and included both 
isovector and isoscalar pairing treated in the LN approximation. The deformation space 
covered quadrupole P2, 7 and hexadecapole P4 shapes. The liquid-drop formula of [^] was 



used to calculate the macroscopic part of the total energy. Apart of including t=0 pairing, 
the details of the method follow our standard implementation and we refere the reader to 



Ref. [0 for further details. 



The calculations have been done for all 0-0 N=Z nuclei from ^^F to ^^Rb, see table |. 
To account for the T = states the 2QP surfaces aJa^lO) have been created where a\(^2) 
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FIG. 4: Pair transfer amplitude (|-Pio|) calculated in the BCS approximation for a model involving 
t=l correlations only. The results are shown as a function of the contribution of Apn to the total 
gap, see text for details. 

denote the lowest qp excitations of the same signature. We use the notion of signature rather 
than time-reversal since the calculations have been done at a small fixed spatial rotational 
frequency. By blocking aJoglO) we assure that all pairing fields are truely blocked. Indeed, 
blocking the lowest qp states of different signatures Q;Ja||0) results in the blocking of isovec- 
tor pp and nn pairing but not isoscalar pn pairing. It is assumed that the valence blocked 
quasiparticles [proton-neutron] may interact. This interaction is added as a perturbation 
assuming a simple [isoscalar] delta force g{A)effS{ri — r2) as a residual interaction. For fur- 
ther simplification we assume either spherical |(n/jm)^; /) or Nilsson asymptotic KA^n^Af^)^) 
limits for the two-body wave function with / = lexp or lexp = '^K, respectively. An effective 
strength of the residual force g{A)eff = geff/VA was assumed. This gives rise to a ~ 1/A 
dependence for the matrix element in accordance to the standard mass dependence of the 
valence pn interaction in liquid-drop models [^ . 



The T = 1 states were calculated in a slightly simplified manner. Self-consistent TRS 
calculations have been done only for the OQP state [false vacuum] at C(j^=0. Then, the 
isocranking calculations were performed at the deformation minimum of the TRS calcula- 
tions. 

The TRS calculations have been performed for a slightly reduced strength of the pairing 
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TABLE I: The lowest T = and T = 1 states in odd-odd N=Z nuclei. Data are taken from: ^^F — 
^^Cu Ref. |§; ^'^Ga Ref. |§; ^^As Ref. 1^; ™5r Ref. H; '^^Kr Ref. |3|. Spin values in square 
parantheses are uncertain. 



force as compared to the estimate of ||14[. The TRS and (5-force corrected TRS results 
are presented in fig. The strength of the 5-force g^ff ~ 650 Me V was estimated from a 
least square fit to the data. The size of the matrix element is more or less consistent with 
average liquid-drop estimate of ~ 20/^4 Me V of the pn effect for valence particles in 0-0 nuclei 
although one one would expect an enhancement in N=Z nuclei purely due to geometrical 
reason (congruence effect). However, the overall quality of the fit does not change very much 
even if we double Qeff and therefore is not very conclusive. 
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FIG. 5: The difference AEt = Et=o — Et=i calculated using the TRS model [triangles]. Open 
circles and diamonds mark the calculations including the correction due to the interaction of the 
valence particles evaluated assuming Nilsson asymptotic and spherical wave functions, respectively 
[see text] . Solid circles mark the available data which are collected in Table |. 



The new elements in these calculations with respect to the one presented in Ref. |]T5 
can be summarized as follows: (z) the equilibrium deformation is calculated from the TRS 
minimum (ii) two quasiparticles of the same type (i.e. not in time reversed states) have 
been blocked and the residual pn interaction between the valence pair has been added to 
the energy difference to reduce the coherence of the t=0 pair field. These new calculations 
describe well the global decrease of AE{A) including the sign inversion in the /7/2 shell, 
see Fig ^. However, the details are far from being satisfactory reproduced. More extended 
studies of the wave function of the valence nucleons is not expected to improve the situation 
since the TRS ground states are essentially spherical for light nuclei. Deformation sets in 
only for heavy nuclei where spherical and asymptotic matrix elements of the 5-force are 
anyhow very similar. This seems to indicate that the t=0 pairing is too strong and needs 
further reduction. 
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VIII. INFLUENCE OF ISOVECTOR PARTICLE-HOLE FIELDS ON THE MO- 
MENT OF INERTIA IN ISOSPACE 



Phenomenological potentials [like the Woods-Saxon potential used here] depend only on 
the third component of the isospin / = {N — Z) / A. It means that essentially no modification 
of the mean-field is done as a function of excitation energy in a given nucleus. On the 
contrary, such modifications are automatically included in self-consistent approaches through 
the changes of the isovector densities. 

The presence of a repulsive isovector mean-field provides an alternative [or additional] 
mechanism to lower the Mol in isospace. Let us illustrate it by using the simple iso-cranked 
single-particle model discussed in detail in Ref. |]14|, |1^ and including an additionally 
repulsive two-body interaction ^ni ■ i analyzed by Neergard [p^ : 



H'^ = h,p- Lji+l-Kt-i (73) 



Linearization of the Hamiltonian ([73|) and simoultanous assumption of one-dimensional rota- 
tion [say around x-axis, then (t^,) = (t^) = 0] leads effectively to a mean-field one-dimensional 
cranking Hamiltonian: 

with an effective isospin dependent cranking frequency. The role of the isovector field is 
depicted schematically in fig. ^ where for simplicity the equidistant single particle (s.p.) 
spectrum Cj = ibt for hsp = Yliocc^i assumed. It is clearly seen from the figure that 
isovector field simply shifts crossing frequencies from 5e, 35e, 55e, ... to 5e+K, 3((5e+K), 5(5e+ 
k), . . . These shifts are marked by dotted lines in fig. |^ since they are in fact instantenous 
due to non-collective iso-cranking. The energy dependence E{T) as a function of isospin 
[T = = (4)] reads: 

E{T) = ]^{5e + 2K)T^ (75) 



i.e. is analogical to the one derived in |T|] but with a reduced Mol-i. Adding standard 
t = 1, = ±1 pairing with A„ = Ap deformes the system [see Sec. |T1[ and smoothes out the 
single-particle, step-like alignment process [collective iso-rotation] but does not essentially 
change the Mol-i as already discussed in |T^, |T3|, |TB| [see also Sec. |111B |. The collectivity 
effect introduced by the t = l,tz = ±1 pairing correlations allows for incorporating isospin 
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fluctuations in a static way through standard cranking condition 

E=^{6e + 2k)T^ ^ 1 (5e + 2k)T{T + 1) (76) 

and thus restores the hnear [Wigner] term. This regime of the model with fixed A, and 
= can be called the vertical excitation regime. Changing the iso-cranking generator 
from to tz and limiting ourselves to total particle number conservation A only brings our 
model to the regime of horizontal excitations describing ground states of neighbouring nuclei 
of the same A. The mathematics used to solve both models is identical. Note, however, that 
the physics interpretation changes. For example, the physical restrictions for allowed s.p. [or 
quasiparticle] excitations which were related to iso-signature conservation and time reversal 



symmetries |jTj, |T5|, now can be interpreted in terms of neutron and proton number 
parities. The cranking frequency measures the difference between neutron and proton Fermi 
energies. 

In the derivation of eqs. ( [75D the isovector Hartree potential Vt = K,{i)i was used rather 
than the two-body interaction. This seems to be consistent with standard potential-like 
treatment of the isovector terms in phenomenological nuclear potentials. In such a case 
use of the cranking condition (^) also makes sense. Since for large isospins pn pairing is 
irrelevant, the inertia parameter [the symmetry energy strength agym] is fully determined by 
the mean-level density and k. The data does not give any signature of enhancement of agym 
in N=Z nuclei. Still isoscalar pra-pairing effects may be related to the enhancement of the 
linear term to \asymT{T + x) with x > 1. 

The symmetry energy strength, and indirectly the value of k, can be conveniently fitted 
using the so-called double difference Vpn ~ dmz formula [^]. The results of a local fit to 
N Z nuclei involving 2^ > 10, 1 < < 3 nuclei [except o-o Tz=l nuclei] are shown 
in table |I|. This fit assumes Esym oc (A^ — Z)^ since the linear term vanishes because 
of the second derivative. However, the strength of the linear term [or the Wigner energy 
E^ig oc — Z|] can be fitted separately using the prescription of Ref. [Q which gives 



information about the local enhancement factor x in the Egym ~ T{T + x) formula. 

The obtained strength of the symmetry energy lies in between large scale fits of Esym ~ 
type [Q and with those assuming Egym. ~ T{T + 1). The latest gives: 

E..^ . (l34.4 - Iil±i) MeV . \ml<^ MeV (77) 
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6g 5e+g \X 36g 3(de+K)\ 



FIG. 6: Single particle routhians representative for a model described by the Hamiltonian (74). 
Filled circles mark occupied states. Arrows indicate shifts in crossing frequencies due to the isospin 
dependence introduced via the i/tt • t interaction. 



0.95 39 0.196 31 0.106 1.26 

1/2 8 0.239 6 0.153 1.33 

2/3 14 0.213 11 0.125 1.27 

1 47 0.196 38 0.107 1.24 

TABLE II: Results of the least-square fit for the Wigner energy strength E^ig = aw\^ — Zl/A" 
and for the symmetry energy strength Egy^ = asym{N — Z)"^ /A°' . The minimum of the mean 
standard deviation is reached at a ~ 0.95 for both Wigner energy and symmetry energy fits. 
Results of one-dimensional fits at fixed a = 1/2; 2/3; 1 are given for comparison. 
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for A ~ 50 [Q. The data show clear enhancement of the linear term with x ~ 1.25 [last 



column] which is consistent with the early findings of Janecke Indeed, it clearly leaves 
room for isoscalar pairing since isospin fluctuations in the static theory limit gives a; = 1 by 
definition. 

We have demonstrated recently that the schematic interaction ( |T5D captures many 
features of realistic effective isovector interactions. In particular, in the Hartree-Fock limit, 
it gives rise to a linear term in the symmetry energy: 

^ \kT. (78) 

This term is only due to the isovector part of mean HF potential and therefore represents 
only a fraction of the linear term. The mean-level density related contribution 

^.(AP) (79) 

seems to go beyond the mean-field approximation and its microscopic evaluation requires 
RPA calculations [|2§ . 



Standard pp/nn pairing leads to a strong quenching of the linear term (|7^), as compared 
to its sp estimate Since for a schematic force k is fixed the suppression of this term 



is entirely due to enhancement of the T/ (Ai^) ratio. This mechanism may also quench the 
linear term (|7^). Therefore, from microscopic point of view, one would expect x < 1 i.e. 
below the static estimate. For example, Neergard [^] gives x ~ 0.8 for T = 2, A = 48 nuclei. 



Let us point out, however, that his estimate is based on a mean- level splitting deduced from 



Fermi gas model which is unrealistically large as shown in |^ . 

Moreover, Neergard's model does not take into account t=0 pairing correlations. These 
correlations, even by entering at the dynamical RPA level, are expected to reduce isospin 
fluctuations as they in fact do in the static BCS or LN theory, see fig. 0. In conlusion, 
as expected, the presence of an isovector particle-hole field will definitely reduce isoscalar 
pairing effects but it does not rule out their existence in iV ~ Z nuclei. 

IX. SUMMARY AND CONCLUSIONS 

We present analytical non-selfconsistent solutions for a model including schematic isovec- 
tor and isoscalar pairing and discuss in detail the response of the t=l and t=0 pair-fields 
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FIG. 7: Dispersion in isospin (At^) calculated for the LN and BCS models as a function of 
the ratio of the strengths of isoscalar to isovector pairing forces. Note that t=0 pairing strongly 
quenches (At^). In particular, in the limit of pure isoscalar pairing (BCS) (At^)=0. 

to rotations in isospace in N=Z nuclei. We are particularly interested in the relations of 
the gauge angles of the t=l pair-gaps and the position of cranking axis in isospace. These 
relations decide upon the character of rotation in isospace. In particular, it is shown that 
within the standard model including nn and pp pair correlations and under the assumption 
of eq.(0) no tilted solutions are possible. Isorotation is either of collective type A ± or 



non-collective type A || uj. Since isorotation in N=Z is planar, the pn t=l pair-field always 
induces collectivity. Again it is shown that the most general solutions are obtained for pure 
collective cases when the t=l pair-field [Aj^, Ap„] is perpendicular to the isocranking axis 
i.e. when A^ -L OJ- Other solutions might be possible only for particular values of the gaps. 

We also demonstrate that, within the simple pair- models, mixing of t=l and t=0 corre- 
lations is essentially forbidden or more precisely restricted to a very special combination of 
pair gaps. This is due to the different phase structure of eigenvectors among time-reversed 
states, compare eq. ([T0|) and eq. (^3|), or alternatively, due to the different transformation 
properties of t=l and t=0 gaps under time-reversal. 

Numerical calculations show that this mixing is allowed within the HFB approximation 
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but only for more elaborate forces [|5|, |TT|]. It is also possible for schematic forces, within 
the LN approximation 0. In the latter approximation, however, the pn t=l pair gap 
vanishes which makes it unsuitable to calculate the T=l pair transfer between the T=l o-o 
ground state and T=0 e-e vacuum. Let us point out very clearly that, since T=l o-o ground 
state and T=0 e-e vacuum are according to our interpretation both OQP states, the t=l pair 
transfer can in principle become enhanced through t=l pn correlations as it is demonstrated 



within the BCS approximation in Sec. 0. On the other hand, the different fundamendal 
structure of the T = vaccua in e-e and o-o nuclei results in a quenching of the T=0 pair 
transfer even in the presence of strong t=0 pairing. 

The energy difference of AEt = Et=q — Et=i can be reproduced in a schematic manner by 
our extended TRS calculations, including a residual pn-force of 5-type. Some discrepancies 
remain, most likely related to the fact that the t=0 pairing is too strong in our calculations 
due to the lack of particle-hole isovector interaction. Inclusion of such an interaction will 
result in a weakening of the isoscalar pairing. However, the enhancement in binding energy in 
N=Z as compared to N-Z=2 nuclei seems to leave quite some room for t=0 pair correlations 
as discussed in Sec. |VII4 Further work along this line is in progress. 

This work was supported by the Goran Gustafsson Fundation, the Swedish Institute, the 
Polish Committee for Scientific Research (KBN) under Contract No. 5 P03B 014 21, and 
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of Tennessee), and DE-AC05-00OR22725 (UT-Battelle). 
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